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ABSTRACT 

We present the results of X-ray observations of four bright transients sources detected in the July 2004 XMM- 
Newton observations of the central bulge of M31. Two X-ray sources, XMMU J004315.5-1-412440 and XMMU 
J004144.7-1-411110, were discovered for the first time. Two other sources, CXOM31 J004309.9H-412332 and 
CXOM31 J004241.8-1-411635, were previously detected by Chandra. The properties of the sources suggest 
their identification with accreting binary systems in M31. The X-ray spectra and variability of two sources, 
XMMU J004144.7-H411110 and CXOM31 J004241.8-H411635, are similar to that of the Galactic black hole 
transients, making them a good black hole candidates. The X-ray source XMMU J004315.5-1-412440 demon¬ 
strates a dramatic decline of the X-ray flux on a time scale of three days, and a remarkable flaring behavior on 
a short time scales. The X-ray data on XMMU J004315.5-1-412440 and CXOM31 J004309.9-1-412332 suggest 
that they can be either black hole or neutron star systems. Combining the results of 2000-2004 XMM observa¬ 
tions of M31, we estimate a total rate of the bright transient outbursts in the central region of M31 to be 6 - 12 
yr ', in agreement with previous studies. 

Subject headings: galaxies; individual (M31) — X-rays; binaries — X-rays; stars 


1. INTRODUCTION 

Bright X-ray transient sources provide a unique opportunity 
to study the properties of the accretion onto stellar-mass com¬ 
pact objects. The observations of the Galactic X-ray novae 
(XRNe) show that they can b e either high-mass (HMXB ) and 
low-mass (LMXB) binaries dTanaka & Shibazaki 19961) . but 
most of confirmed black-hole binaries are low-mass systems 
(McClintock & Remillard 2005, and references therein). 

Until recently, the detailed study of X-ray transient sources 
has been mostly limited to our Galaxy. The advent of a new 
generation of X-ray telescopes {Chandra and XMM-Newton) 
has allowed the study of the spectral and temporal properties 
of XRNe located in the nearby galaxies and compare them to 
their Galactic counterparts. 

The relative proximity and favorable orientation of M31 
make it a prime target for the study of an extragalactic XRNe 
population. Recent Chandra and XMM-Newton observations 
of M31 have led to the discovery of several dozen bright tran¬ 
sient X-ray sources with luminosities between ~ 10^® and 
~ 10^* ergs s“' (Garcia et al. 2000; Trudolyubov et al. 2001; 
Kong et al. 2002; Williams et al. 2004; Williams et al. 2005a, 
b, c). The follow-up observations of X-ray transients with 
HST ensured the identification and study of optical counter¬ 
parts to some of these systems (Williams et al. 2004; Williams 
et al. 2005a, b). 

Here we present the results of spectral and timing analy¬ 
sis of four bright transient sources detected in the 2004 July 
XMM-Newton observations. 

2. OBSERVATIONS AND DATA ANALYSIS 

The central region of M31 was observed with 
XMM-Newton on four occasions during July 2004 
llBarnard et al. 200^ (Eig. Q. In the following analysis 

* Institute of Geophysics and Planetary Physics, University of California, 
Riverside, CA 92521 
Electronic address: sergeyt@ucr.edu 

^ Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 
84/32, Moscow, 117810 Russia 

^ Los Alamos National Laboratory, Los Alamos, NM 87545 


we use the data from three European Photon Imag¬ 
ing Ca mera (EPIC) instru ments; two EPIC MOS de- 
tectors j Turner e t al. 200lb and the EPIC-pn detector 
llStrueder eTaT^OOll) . In all observations EPIC instruments 
were operated in the full window mode (30' EOV) with the 
medium optical blocking filter. 

We reduced EPIC data with the latest version of XMM- 
Newton Science Analysis System (SAS v 6.5.0)^^. Each of the 
original event files were screened for periods of high back¬ 
ground. The remaining exposure times for each observation 
are listed in Tabled The 2004, July 17 observation (Ohs. #2 
in Table [3 is affected by high background, so we excluded 
it from our spectral and timing analysis, and used it for total 
flux estimates only. 

We generated EPIC-pn and MOS images of the central re¬ 
gion of M31 (Pig. in the 0.3 - 7.0 keV energy band, 
and used the SAS standard maximum likelihood (ML) source 
detection script edetectjchain to detect and localize point 
sources. We used bright X- ray sources with k nown optical 
counterp ^ts from USNO- B JMonet et al. 2003h and 2MASS 
catalogs JCutri et al. 2003h to correct EPIC image astrometry. 
After applying the astrometric correction, we estimate resid¬ 
ual systematic error in the source positions to be of the order 
0.5 - 1". To identify transient sources, the resulting list of 
the XMM sources was compared to the existing catalogs of 
M31 X-ray sources (Trinchieri & Pabbiano 1991; Primini et 
al. 1993; Kong et al. 2002; Williams et al. 2004; Pietsch et al. 
2005) and transient source lists from the Chandra monitoring 
campaign (e.g. Williams et al. 2005a, b, c). 

To generate lightcurves and spectra of X-ray sources, we 
used elliptical extraction regions with semi-axes size of ~ 
15 - 50" (depending on the distance of the source from the 
telescope axis) and subtracted as background the spectrum of 
adjacent source-free regions, with subsequent normalization 
by ratio of the detector areas. Por spectral analysis, we used 
data in the 0.3 - 7 keV energy band. All fluxes and luminosi¬ 
ties derived from spectral analysis apply to this band. We used 
spectral response files generated by XMM SAS tasks. Spectra 
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were grouped to contain a minimum of 20 counts per spec¬ 
tral bin in order to allow statistics and fit t o analytic mod - 
els using the XSPEC v.ll^ fitting package JArnaud 19961) . 
EPIC-pn, MOS1 and MOS2 data were fitted simultaneously, 
but with normalizations varying independently. The energy 
spectra of the sources were fitted by two standard X-r ay bi¬ 
nary spectral models jMcClintock & Remillard 2005|) : ab¬ 
sorbed simple power law or multicolor disk blackbody models 
(DISKBB). To estimate upper limits on the quiescent source 
luminosities, the EPIC count rates were converted into energy 
fluxes in the 0.3 - 7 keV energy band using Web PIMMS®, 
assuming standard parameters: an absorbed simple power 
law m odel with Ah - 7 x 10^^ cm“^ (Galactic foreground 
value (iDi ckev & Lockman 19901) 1 and photon index a = 1.7 
dShirev et al. 200ll) . We used standard XANADU/XRONOS 
v.5^ tasks to perform analysis of the timing properties of the 
transient X-ray sources. 

In the following analysis we assume M31 distance of 760 
kpc (van den Bergh 2000). All parameter errors quoted are 
68% (Icr) conhdence limits. 

3. RESULTS AND DISCUSSION 
3.1. XMMUJ004144.7+411 no 

A new X-ray source XMMU J004144.7H-411110 has been 
discovered in the data of 2004 July 16 XMM-Newton observa¬ 
tions of M31. In addition, our analysis of the archival Chan¬ 
dra data (2004, July 17 observation #4719) revealed the pres¬ 
ence of a bright X-ray source at the position consistent with 
XMMU J004144.7H-411110. Combining the data of XMM- 
Newton and Chandra observations, we measure the position 
of XMMU J004144.7H-411110 to be a = 00''41'"44.70^ d = 
41°11'lO" (J2000 equinox) with an uncertainty of ~ 1". The 
search for the optical coun terparts using the im ages from Lo¬ 
cal Group Survey (LGS) (iMassev et al 200ll) did not yield 
any object brighter than ~ 21 within the error circle of 
XMMU J004144.7H-411110. Using the data of archival XMM 
observations of the central region of M31, we estimate the 
upper limit {2cr) on the source quiescent luminosity to be 
~ 2x10^^ ergs s“' in the 0.3-7 keV energy band, >100 times 
lower than maximum measured outburst luminosity (Table|2}. 

The energy spectra of XMMU J004144.7H-411110 are soft, 
and can be well htted by absorbed DISKBB model with color 
temperatures ~ 0.6 - 0.8 keV or absorbed power law model 
with photon index of ~ 2.8 - 3.3 (Tabled Eig. |2j. The 
corresponding estimated luminosities of the source have been 
found to be in the range of ~ (2.3-3.0)xl0^^ ergs s“'. Eortwo 
observations (1 and 3), the energy spectrum shows clear signs 
of high energy cut-off: the DISKBB model approximates it 
better than a simple power law, as indicated by fit statistics 
(TableEJ. 

The X-ray spectrum, transient behavior and extreme 
faintness of the optical counterpart indicate that XMMU 
J004144.7-I-411110 is not a Galactic foreground object, and 
probably is an accreting binary system in M31. The spec¬ 
tral model hts require absorbing columns well in excess the 
Galactic foreground value of 7 x 10^** cm“^ (Table |2j, this 
could be cons istent with the source located insid e or behind 
the M31 disk dTrudolyubov & Priedhorsky 2004|) . 

The observed spectrum and luminosity of XMMU 
J004144.7-I-411110 bear clear resemblance to the Galactic 

^ http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 
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black-hole transients in the /u'g/j/“thermal-dominant” state 
during the flux decline that precedes the transition to the 
lowlhard state (Tomsick & Kaaret 2000; McClintock & 
Remillard 2005). It should be also noted, that the 0.3 - 7 keV 
spectrum of the source is signihcantly softer than observed 
in the neutron star systems at similar lu minosity levels both 
in the G alaxy ( Christian &_Swank 1997 |) an d in M31 globular 
clusters dTrudolyubov & Priedhorsky 2004ll . 

3.2. XMMU J004315.5+412440 

2004 July XMM observations revealed another previously 
undetected X-ray source, located at a = 00^'43'”15.5U, d = 
41°24'40" + 1.5". The inspection of the LGS images showed 
no optical sources brighter than my ~ 21 in the XMM error 
circle of XMMU J004315.5-1-412440. Using the data of pre¬ 
vious XMM observations, we estimated quiescent luminosity 
of the source to be < 10^^ ergs s“' in the 0.3 - 7 keV energy 
band. 

The X-ray source XMMU J004315.5-H412440 demon¬ 
strates a remarkable variability on time scales ranging from 
minutes to several days (Eig. |3l. The source flux in the 
0.3 - 7 keV band changed dramatically in the course of four 
2004 July XMM observations, dropping from ~ 10“^^ erg s“* 
on July 16 to ~ 1.5 x 10“'"^ erg s“' in four days (Eig. |3 
upper panel). The source XMMU J004315.5-1-412440 also 
shows a high level of variability during July 16 observation 
(Obs. l)(Fig. El lower panel). In addition to the irregu¬ 
lar flux variations, the source produced an intense flare at 
~ 13700 s (Fig. El lower panel), lasting for ~ 1000 s. The 
time evolution of the source flux during the flare is charac¬ 
terized by fast rise (< 200 s) to a maximum level followed 
by quasi-exponential decay with estimated e-folding time of 
~ 500 + 200 s. The peak intensity of the flare was > 4 times 
higher than the average source intensity, corresponding to a 
luminosity of ~ 3 x 10^^ erg s“' in the 0.3 - 7 keV energy 
band. The estimated total energy emitted during the flare is 
~ 10"^° ergs. Unfortunately, the sensitivity of our observations 
does not allow to make a reliable conclusion on the evolution 
of the source spectrum during this flare. 

The X-ray spectrum of XMMU J004315.5 h- 4 12440 during 
July 16 observation (Obs. 1) is soft, and can be well approxi¬ 
mated by an absorbed power law model with photon index of 
~ 3.8 or by DISKBB model with color temperature of ~ 0.36 
keV (TableEJ- The corresponding absorbed luminosity of the 
source was ~ 7 x 10^® erg s '. The power law model provides 
better fit to the observational data than a DISKBB model, as 
seen from Table El The power law model fit requires a high 
level of low-energy absorption ~ 5x10^' cm“^, while the 
DISKBB model requires an absorbing column consistent with 
Galactic foreground value in the direction of M31 (Table El- 
We did not detect a statistically significant change in the shape 
of the source spectrum during the overall flux decline, as seen 
from power law model fits (Table E}- 

The X-ray properties of XMMU J004315.5H-412440 along 
with the constraints on the optical counterpart support its 
identification as an accreting binary system in M31. A 
combination of strong long and short-term variability of X- 
ray flux makes XMMU J004315.5-1-412440 especially in¬ 
teresting. The profile and luminosity of the X-ray flare 
detected in the July 16 observations are somewhat simi¬ 
lar to the long thermonuclear X-ray bursts detected from 
the Galactic neutron stars, although the spectrum of the 
source is signihcantly softer than that of the burst sources 
l lStrohina'^r^ Bildsten 2005]) . The other possibility is that 
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the flare is caused by spasmodic accretion onto the com- 
pact object simi lar to the neutron star Type II bursts 
(iLewin et al. I995I) . 

3.3. CXOM31 J004241.8+411635 

The X-ray source CXOM31 J004241.8H-411635 was first 
detected in the July 17, 2004 Chandra observation and 
remained detectable i n the 2004 Sept. 2 and Oct. 4 
Chandra observations dWilliams et al. 2005(J) . The source 
has been previously detected in 1979 Einstein observations 
dTrinchieri & Fabbiano 199T1) , which makes it a probable re- 
current transient cand idate with duty cycle of 0.02 - 0.06 
dWilliams et al. 2005d) . In addition, the follow-up observa¬ 
tions of the source with HST resulted in a determination that 
the optical counterpart to CXO M31 J004241.8-H41163 5 is a 
low-mass binary system in M31 dWilliams et al. 2005(Jl . 

The X-ray spectra of CXOM31 J004241.8H-411635 mea¬ 
sured with XMM/EPIC are relatively soft (Fig. |2j, and can be 
approximated by absorbed power law with photon index of 
~ 2.3 - 2.5, or by DISKBB model with characteristic temper¬ 
ature of ~ 0.86 - 0.99 keV (TableE). Due to the curvature in 
the source spectrum around 1 keV, the DISKBB model pro¬ 
vides significantly better description to the data than a power 
law. The power law model fits result in a relatively high val¬ 
ues of the absorbing column (A^h ~ 2 x 10^' cm“^), while 
in the DISKBB approximation, the measured absorption was 
consistent with Galactic foreground absorption toward M31. 
The X-ray spectra can be also fitted by combination of the 
DISKBB and power law models with kTjn ~ 0.9 - 1.1 keV 
and a photon index of ~ 4-5, improving the overall quality 
of the fit with respect to one-component model approxima¬ 
tion. 

The estimated luminosity of the source changes from 5 x 
10^’ to ~ 6.5 X 10^^ ergs s“' in the course of XMM observa¬ 
tions. There is a clear correlation between the spectral temper¬ 
ature derived in the DISKBB spectral fits and the X-ray flux 
(Table 13, sim ilar to the results of th e Chandra observations 
of the source dWilliams et al. 2005(J) . At the same time, the 
characteristic emitting radius /?,„ sjcos i remains essentially 
constant, despite a significant change in the X-ray flux (Table 
13, the effect also observed in the Galactic black hole candi¬ 
dates in the high spectral state l | Tanaka & Lewin 1995b . 

The X-ray properties of CXOM31 J004241.8H-411635 ob¬ 
served with XMM have been found to be remarkably similar to 
the properties of the Galactic black hole candidates in the hig h 
or “thermal-dominant” state dMcClintock & Remillard 2005L 
in ge neral agreement with t he results of Chandra observa¬ 
tions dWilliams et al. 2005<jl . which makes it a good black 
hole candidate. 

3.4. CXOM31 J004309.9+412332 

The X-ray source CXOM31 J004241.8H-411635 was 
first detected in the 2004 May 23 Chandra observation 
dWilliams et al. 2005M) . Regular Chandra monitoring of the 
source revealed a complex X-ray light curve with at least 
two peaks separated by ~ 130 days and reaching the max¬ 
imum 0.3 - 7 keV luminosity of ~ 10^^ ergs s“^ Based 
on the results of Chandra and HST observations, CXOM31 
J004309.9-H41 2332 has been classifie d as a low-mass X-ray 
binary system dWilliams et al. 2005 Id) . 

The XMM-Newton spectra of CXOM31 J004241 .8 h- 41 1635 
during 2004 July XMM observations can be approximated by 
steep absorbed power law with photon index of ~ 2.8 - 4.0 


or by DISKBB model with kr,„ ~ 0.33 - 0.36 keV (Table 
13. The corresponding absorbed luminosity of the source 
is in the range of ~ (2 - 3) x 10^® ergs s“’. The spec¬ 
tra of CXOM31 J004241.8-1-411635 measured with XMM- 
Newton, correspond to the decline of the first outburst, 
and appear to be significantly softer than measured during 
2004, May 23 and Sep. 2 C handra observations (c oinci- 
dent with two outburst peaks) dWilliams et al. 2005^11 . The 
anti-correlation of the spectral hardness and X-ray flux of 
CXOM31 J004241.8-H411635 could be similar to that ob¬ 
served in some Galactic low-mass X-ray binaries (both neu¬ 
tron star and black h ole binary systems) and in globuly clus¬ 
ter sources in M31 dTriidolviihov & Priedhorskv 2004[) . The 
available data on CXOM31 1004241.8-1-411635 suggest that 
it could be either a black hole or a neutron star, not allowing 
to make a definitive conclusion on the nature of the compact 
object in this system. 

4. CONCLUSIONS 

Using the data of 2004 July XMM/EPIC observations of 
M31, we study the X-ray properties of four bright transient 
sources. Two X-ray sources, XMMU J004315.5-1-412440 
and XMMU J004144.7-1-411110, were discovered for the 
first time. Two other sources, CXOM31 J004309.9-1-412332 
and CXOM31 J004241.8-1-411635, were previously detected 
by Chandra. The properties of the sources along with 
the information on optical counterparts suggest their iden¬ 
tification with accreting binary systems belonging to M31. 
The luminosities, energy spectra and variability of the 
two sources, XMMU J004144.7-I-411110 and CXOM31 
J004241.8-1-411635, are reminiscent of the Galactic black 
hole X-ray novae, which makes them a good black hole candi¬ 
dates. The X-ray source XMMU J004315.5-1-412440 demon¬ 
strates a dramatic decline of the X-ray flux on a time scale of 
three days, and a remarkable flaring behavior during the first 
XMM observation. The available data on this source and the 
other transient source CXOM31 J004309.9H-412332 is con¬ 
sistent with either black hole or a neutron star interpretation. 

The total of 10 transient sources with 0.3 - 7 keV luminosi¬ 
ties higher that 10^® ergs s“^ have been detected in 5 XMM- 
Newton observations of the central part of M31 (Osborne et 
al. 2001; Trudolyubov et al. 2001; Pietsch et al. 2005; this 
Letter). More than a half of these sources (60%) can be classi¬ 
fied as black hole candidates. The remaining sources include 
two supersoft transients with probable white dwarf primaries 
(20%), and two systems, which can contain either black holes 
or neutron stars. 

Given a number of the detected sources, and the cover¬ 
age factor of the XMM observations, one can estimate the 
expected total rate of X-ray transient outbursts in the bulge 
and inner disk of M31. Assuming the average duration of a 
bright phase of the typical transient of ~ 1-2 months, gives 
a total rate of ~ 6 - 12 outbursts per year^, consistent with 
estimates based on earlier XMM-Newton and Chandra results 
(Trudolyubov et al. 2001; Kong et al. 2002; Williams et al. 
2004). 

5. ACKNOWLEDGMENTS 

Support for this work was provided through NASA Grant 
NAG5-12390. XMM-Newton is an ESA Science Mission 
with instruments and contributions directly funded by ESA 

* Note a large uncertainty of this estimate, due to the limited statistics and 
scatter in the observed outburst durations. 






























4 


TABLE 1 

2004 July XMM-Newton observations of the central region of M31. 


Ohs. # 

Date, UT 

UT 

(h:m:s) 

Ohs. ID 

RA (J2000)“ 
(h:m:s) 

Dec (32000)“ 
(d:m:s) 

Exp.(pn)'’ 

(ks) 

#1 

2004 Jul 16 

16:17:05 

0202230201 

00:42:42.12 

41:16:57.1 

16 

#2" 

2004 Jul 17 

12:07:53 

0202230401 

00:42:42.08 

41:16:56.9 

18 

#3 

2004 Jul 18-19 

23:49:30 

0202230401 

00:42:42.26 

41:16:58.2 

13 

#4 

2004 Jul 19 

12:48:23 

0202230501 

00:42:42.23 

41:16:57.6 

6.5 


“ - coordinates of the center of the field of view 


- instrument exposure used in the analysis 

- observation affected by high background 


Member states and the USA (NASA). This research has made Science Archive Research Center Online Service, provided 
use of data obtained through the High Energy Astrophysics by the NASA/Goddard Space Flight Center. 
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TABLE 2 

Model fits to the energy spectra oe transient sources (XMM/EPIC data, 0.3 - 7 keV energy range). 


Obs. 

Model 

Nh 

kT 

Ri„ yjcos P 

Photon 

Flux* 



Instrument 

# 


(xlO^® cm" 

-2) (keV) 

(km) 

Index 


(d.o.f) 



XMMU J004144.7+411110 

1 

PL 

51±5 



3.16 ±0.15 

4.38 ±0.16 

56.6(53) 

30.3 pn+Ml+M2 


DISKBB 

19 ±3 

0.68 ± 0.04 

30 ±4 


4.32 ±0.16 

48.8(53) 

29.9 


3 

PL 

40y 



9 79+0.14 
^•'^-0,12 

3.59 ±0.16 

64.8(53) 

24.8 pn+Ml+M2 


DISKBB 

14 ± 3 

n '70+0.06 

20!' 


3.36 ±0.15 

67.6(53) 

23.2 


4 

PL 

55 ±7 



Q Q 1 +0.12 

J.J1_o22 

3.89 ± 0.20 

54.9(44) 

26.9 

pn+M2 


DISKBB 

22 ±4 

0.60 ± 0.06 

3T3, 


3.66 ±0.19 

51.0(44) 

25.3 


XMMU J004315.5+412440 

1 

PL 

29 ±5 



Q 01 +0.44 
■^•°^-0.32 

1.07 ±0.06 

33.9(22) 

7.39 

pn 


DISKBB 

'7+4 

0 

^••^^-0.04 

51 ± 14 

0.98 ± 0.06 

38.7(22) 

6.77 


3 

PL 




-1.67 

0.20 ± 0.03 

4.1(5) 

1.38 

pn 

CXOM31 J004241.8+411635 

1 

PL 

23 ±2 
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i. 1. — The central region of M31 as it appears in XMM-Newton EPIC-MOS 2004 July 16 observation. The transient sources are marked with arrows. 
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Fig. 2.— Count spectra of bright transient sources. EPIC-pn data, 0.3 - 7 keV energy range. The analytical model fits are shown with thick 
histograms. Upper left: Transient source XMMU J004144.7+411110 (absorbed DISKBB model). Upper right: XMMU J004315.5+412440 
(absorbed power law model). An absorbed power law model fit to the data of July 18 observation is shown with dotted histogram. Lower 
left: CXOM31 J004241.8+411635 (absorbed DISKBB+power law model). Lower right: CXOM31 J004309.9+412332 (absorbed power law 
model). 
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Fig. 3.— a: The long-term X-ray light curve of the transient source XMMU J004315.5+412440 during four 2004 July XMM-Newton observa¬ 
tions, obtained from combined data of EPIC-pn, MOSl and MOS2 cameras in the 0.3 - 7.0 keV energy range, b: Background and vignetting 
corrected X-ray light curve of XMMU J004315.5+412440 during 2004 July 16 XMM-Newton observation (Obs. #1), obtained from data of 
EPIC-pn in the 0.3 - 3.0 keV energy range, with a 300 s time resolution. 





